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Priifungsantrag gem. § 44 PatG ist gestellt 

© Computertomograph mit reduzierter Dosisbelastung bzw. reduziertem Bildpunktrauschen 

© Die vorliegende Erfindung schafft einen Computerto- 
mograph (Ein- oder Mehrzeiler mit oder ohne Spiralabta- 
stung), bei dem das Rauschniveau der interpolierten Pro- 
jektionen einen bestimmten Schwellwert nicht ubersteigt 
und der bei gleichbleibender Bildqualitat eine signifikante 
Reduzierung der Patientendosis erlaubt oder dement- 
sprechend bei gleichbleibender Patientendosis eine signi- 
fikant hohere Bildqualitat erzielt. 

Die zu rekonstruierenden planaren Projektionsdaten wer- 
den nach folgender Formel berechnet: 

die Filterung findet also entweder vor oder nach Umrech- 
nung auf Parallelgeometrie statt und minimiert somit den 
™ glattungsbedingten Verlust an Ortsauflosung. 
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Beschreibung 

Die modeme Computertomographie (CT) beinhaltet eine Vielzahl von MeBmodalitaten. Einerseits existiert die kon- 
ventionelle CT, also die Aufnahme einzelner Schichten. Diese wurde groBtenteils durch die Volumenaufnahmetechnik 
5 Spiral-CT ersetzt. Andererseits werden neuerdings CT-Gerate gefertigt, die statt nur einer Detektorzeile mehrere Zeilen 
(N Zeilen) besitzen. Mit diesen Mehrzeilencomputertomographen sind sowohl Schichtaufnahmen von nunmehr N 
Schichten gleichzeitig, als auch Spiralaufnahrnen moglich. 

All diesen Modaliuilen geiiieinsain isu daB eine Konvenlion ("z-Inlerpolation") X beslelu, die gcniessenen Dalen pflj, 
a, n) ((3 ist der Winkel im Facher, a der Projektionswinkel also die Winkelsiellung der Rohre bzw. Winkelstellung von 

10 Rohre und Detektorsystem und n zahlt die Nummer der Detektorzeile oder die Nummer der Schichtaufnahme) auf einen 
planaren transaxialen an der z-Position z R giiltigen Datensatz px(P, 0t, z R ) (Schichtdatensatz) umzurechnen. Dieser 
Schichtdatensatz wird dann mit einem Rekonstruktionsverfahren fiir planare Daten (ublicherweise gefilterte Riickprojek- 
tion oder Fourierrekonslruktion ! rekonstruiert, urn das CT-Bild /u erhalten. 

Bei konventionellen Aufnahmen (Schichtaufnahmen). die aus Einzelumlaufen an Positionen /', z". /"', . . . entstehen, 

15 wird iibhcherweise aber nicht notwendigerweise die Rekonstruktion nur an genau diesen Positionen zugelassen (d. h. z R 
e {z', z", z'", ...}). Bei Spiralaufnahmen ist die Rekonstruktionsposition z R im allgemeinen frei und retrospektive wahl- 
bar. Die am weitesten verbreitete Konvention zur z-Interpolation bei Einzeilen-Spiral-CT ist der Algorithmus X = 
180°LI, cine lincarc Interpolation zwischcn in z-Richtung gcmcsscncn Datcnpunktcn [Willi A. Kalcndcr Wolfgang 
Seissler, Ernst Klotz und Peter Vbck, "Spiral volumetric ( T with single-breath-hold technique, continuous transport, and 

20 continuous scanner rotation", Radiology 176 (1), S. 181-183, Juli 1990], [Arkadiusz Polacin, Willi A. Kalender und Guy 
Marchal, "Evaluation of section sensitivity profiles and image noise in spiral CT", Radiology 185 (1), S. 29-35, Oktober 
1992]. 

Bei alien CT-Bildern gilt, daB die Bildqualitat in bezug auf Rauschen und Niedrigkontrasterkennbarkeit monoton mit 
der Paticntcndosis stcigt. Das hciBt, daB der Anwcndcr den Rohrcnstrom und somit die Dosisbclastung des Paticntcn so- 

25 weit erhoht, bis die Bildqualitat seinem (subjektiven) Empfinden nach "gut" ist. Dosisbetrachtungen haben insbesondere 
in europaischen Fandcrn einen liohen Slellenweil liikI soniil werden \ iele Versuehe unlernonnnen. die Paticntcndosis bei 
moglichst gleichbleibender Bildqualitat zu reduzieren. 

Eine einfache Moglichkeit zur Reduzierung des Bildpunktrauschens bietet sich durch die Wahl eines glattenden Re- 
konslruklionslillers. Die Rekonslruklionslilter konnen anjedem kommerziell erwerblichen CT-Gerat in gewissen Gren- 

30 zen frei gewahlt werden und somit kann jeder Anwender das Rauschniveau im Bild verringern, ohne die Paticntcndosis 
zu erhohen. Da mit dieser Methode der gesamte Datensatz geglattet wird, geht dies zu angsliiulig mil einer Vcrschlech- 
terung der Ortsauflosung einher. Das Problem, die Patientendosis bei moglichst gleichbleibender Bildqualitat zu reduzie- 
ren, bzw. die Bildqualitat bei gleicher Patientendosis zu erhohen ist damit nichl gelost, denn es ruuB ein KompromiB zwi- 
schen Ortsauflosung und Bildpunktrauschen gemacht werden. 

35 Zur Losung des Problems finden sich in der Literatur Ansatze zur adaptiven Filterung der MeBdaten, d. h. der Daten- 
satz wird nichl global, sondcrn nur lokal geglattet LJiang Hsieh. "Adaptive trimmed mean for computed tomography 
image reconstruction", Proc. of SPTE, 2299, S. 316-324, 1994], [Jiang Hsieh, "Generalized adaptive median filter and 
their application in computed tomography", Proc. of SPIE, 2298, S. 662-672, 1994], [Jiang Hsieh, "Adaptive filtering 
approach to the streaking artifact reduction due to x-ray photon starvation", Radiology 205 (P), S. 391, 1997], [Berkman 

40 Sahiner and Andrew E. Yagle, "Reconstruction from projections under timefrequency constraints", IEEE Transactions 
on Medical Imaging, 14 (2), S. 193-204, Juni 1995]. Ublicherweise werden zur adaptiven Filterung benachbarte Detek- 
torelemente verwendet die Filterung findet also ausschlieBlich in P-Richtung statt. 

Ansatze, die Filterung rohdatenbasiert adaptiv (angepaBt an jeden einzelnen gemessenen Projektionswert) in alien drei 
Dimensionen P-, a- und z-Richlung) durch/.uluhrcn sind bisher nichl bekannl. 

45 Der in Patentanspruch 1 angegebenen Erfindung liegt die Aufgabe zugrunde, fiir einen Computertomographen (Ein- 
oder Mehrzeilersystem mit oder ohne Spiralabtastung) das Bildrauschen und die storenden Rauschstrukturen, die durch 
Korrelation benachbarter Bildpunkte entstehen, weitestgehend zu reduzieren. Damit wird eine Verbesserung der Bild- 
qualitat bei gleichbleibender Patientendosis erreicht bzw. eine Dosisreduzierung bei gleichbleibender Bildqualitat mog- 
lich. 

50 Diese Aufgabe ist erfindungsgemaB gelost durch die Merkmale des Palenlanspruchs. 
Die Erfindung ist nachfolgend erlautert. 

Es bezeichne px(^, "t), z R ) die planaren Projektionsdaten (z. B. Schwachungswerte) eines CT-Scans, die an der z-Posi- 
tion z R giiltig sind. Dabei ist ^ der Kanalindex und entspricht in Fachergeometrie dem Winkel im Facher und in Parallel- 
geometrie dem Abstand des Strahls zum Rotationszentrum. -d ist der Projektionsindex und entspricht sowohl in Facher- 

55 als auch in Parallelstrahlgeometrie der Winkelstellung der Rohre bzw. der Winkelstellung von Rohre und Detektor. X sei 
der Name der Konvention zur z-Interpoladon, mil der diese planaren Daten aus der Messung erzeugl wurden. Die Re- 
konstruktionsposition z R sei retrospektiv auf einem der Messung entsprechenden Interv all [z^,,, z max ] frei wahlbar. Dazu 
drei Beispiele. Beispiel 1 : Im Falle der Einzeiler-Spiral-CT ist X = 1 80°LI ein denkbarer z-Interpolationsalgorithmus, der 
hier zur Anwendung kommen kann. Beispiel 2: Bei Mehrzeiler-Spiraldaten kann, analog zu Beispiel 1, eine Standard z- 

60 Interpolation zugrunde liegen. Beispiel 3: Bei konventioneller CT (Schichtaufnahmen ) kann sowohl bei Ein- als auch bei 
Mehrzeilern als Konvention X eine Nachster-Nachbar-Interpolation oder eine lineare Interpolation zweier benachbarter, 
gemessener Schichten verwendet werden. 

Bei dem hier beschriebenen Rekonstruktionsverfahren fiir die Computertomographie werden die planaren Projekti- 
onsdaten p x (£, ft, z R ) vor der Rekonstruktion adaptiv nach folgender Formel gefiltert: 

65 

p AF (^, A z) = ld^dt>'dz'f A ^(^')f Af> (f>-t>')f Az (z-z')p x (^, z'). 
Dabei bezeichnen 



2 



4/8/2010, EAST Version: 2.4.1.1 



DE 198 53 143 A 1 

- pxfe 7 ) die vor Durchfiihrung der adaptiven Filterung (AF) zur Verfiigung stehenden Projektionsdaten in Par- 
allel- oder Facherstrahlgeometrie. Diese errechnen sich aus den MeBdalen durch die Konvention X ^ AF fiir belie- 
bige z. 

- Paf(£> z) die nach Durchfiihrung der adaptiven Filterung zur Verfiigung stehenden Projektionsdaten in Paral- 
lel- oder Facherstrahlgeometrie. Diese Daten konnen dann der Rekonstruktion (z. B. gefilterte Riickprojektion) zu- 
gefiihrt werden. 

- AJj, At*} und A7. die Filterweiten in den drei Koordinatenrichtungen. Die Filterweiten sind Funktionen des aktuell 
adaptiv zu filternden Projektionswerts px(^ ft, z): AE, = A£(p x (cj, z)), Ad = Ad(p x (£, d, z)) und Az = Az(p x (£, d, 
z)). 

- tA^( • ), fAfl( • ) und f^C ■ ) die Filterfunktionen (achsensymmetrisch mit Werten > 0 und Gesamtflache 1 ) fiir die 
Glattung in den jeweiligen Koordinaten. Die Filterweiten A%, Ad und Az stehen jeweils fiir die Halbwertsbreiten 
oder ein anderes charakteristisches BreitemaB der Filterfunktionen. Wird eine oder mehrere der Weiten gleich Null, 
so reduziert sich die Filterfunktion zu einer Diracschen Deltafunktion und es findet in den entsprechenden Koordi- 
nalen keine Filterung statt. 

In einer moglichen Ausfiihrung konnte die Wahl der Filterweiten mit der Formel 

gftr«A»> f a l ]s gPxW) > T > 

sonst 



(A£ + £)(A0 + S)(Az + z) = c< 



(A£ + £):(A,9 + ,9):(Az + z) = £ :& :z 

erfolgen. Dahei hezeichnen 

- cj, ■& und z den Abtastabstand in den drei Koordinatenrichtungen: \ ist die Breite eines Detektorelements, d ent- 
spricht der Detektorintegrationszeit und z der Halbwertsbreite des Schichtempfindlichkeitsprofils. 

- e Px( ^' z) die relative Varianz des Projektionswerts px(%, d, z) als MaB fiir das Rauschniveau des Punkts (cj, ft, z). 

- T 2 den Schwellwert ab dem die Filterweiten A^, At) und A/, einen von Null verschiedenen (positiven) Wert an- 
nehmen. 

- c die Normierungskonstante, die so gewahlt wird, daB bei Unterschreilung des Schwellwerls T die h'illerwcileii 
A^, Ad und Az den Werl Null annehmen. 

Folglich gilt cT 2 = dfz. 

Die Bedingung (A£ + 1) : (Ad + d) : (Az + z) = \ : d : z stellt sicher, daB das Verhaltnis der effektiven Filterweiten (dies 
sind die Filterweiten plus der intrinsischen Mittelung durch das CT-Gerat.) dem geratebedingt.en und somit vom Herstel- 
ler optimierten Verhaltnis aus Detektorbreite, Integrationszeit und z-Ausdehnung entspricht. 

Die Filterfunktion selbst kann z. B. eine der folgenden Darstellungen haben: 

„ _ , _ l fl falls 2W< Ax 

• Rechteck: / (*) = _ J 11 

Ax [o sonst 

_ . . , , , 1 \\-\x\lhx falls|x|<Ax 

• Dreieck: fjjt) = — \ 

Ax 0 sonst 



Diese Darstellungen haben den Vorteil, daB die zugehorigen Integrationen analytisch ausgefflhrt werden konnen und 
die adaptive Filterung sich auf eine gewichtete Summe iiber die gemessenen Daten reduziert. Damil lassen sich die ad- 
aptiven Filter effizient implementieren. 

Patentanspriiche 

1. Computertomograph (Ein- oder Mehrzeiler mit oder ohne Spiralabtastung), bei dem durch geeignete Filterung 
sichcrgcslclll wird, daB das Rauschniveau der Projektionen einen beslimmten Schwellwerl nichl ubersteigt. 

2. Computertomograph nach Anspruch 1, bei dem die Interpolation der Projektionen adaptiv erfolgt (adaptive Fil- 
terung, adaptive z- Interpolation). 

3. Computertomograph nach Anspruch 2, bei dem die Projektionsdaten in Parallel- oder Facherstrahlgeometrie so- 
wohl in Kanalrichtung (c^-Koordinate), in Projektionsrichtung (d-Koordinate) als auch in axialer Richtung (z-Koor- 
dinate) adaptiv gefiltert werden. 

4. Computertomograph nach Anspruch 3, bei dem die Projektionen adaptiv nach folgender Formel gefiltert werden 
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(3D-adaptive Filterung): 

p AF fe 0, Z) = Jd^dd'dz'fA^^^Wz-Z^px^', Z'). 

5. Computertomograph nach Anspruch 4, der bei nahezu gleichbleibender Bildqualitat eine signifikante Reduzie- 
rung der Patientendosis erlaubt, oder dementsprechend bei gleichbleibender Patientendosis eine signifikant hohere 
Bildqualitat. erz.ielt. 
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